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Abstract 

We extend some previous results of our work [Ij on the error of the averag- 
ing method, in the one-frequency case. The new error estimates apply to any 
separating family of seminorms on the space of the actions; they generalize our 
previous estimates in terms of the Euclidean norm. For example, one can use 
the new approach to get separate error estimates for each action coordinate. 
An application to rigid body under damping is presented. In a companion 
paper [2] , the same method will be applied to the motion of a satellite around 
an oblate planet. 
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1 Introduction. 



It is often stated by applied mathematicians that a good theorem on differential 
equations is one outlining a computational method for their solutions; if the method 
is approximate, quantitative error estimates should be provided. 

In the case of ODEs with slow variables ("actions") and fast angular variables, 
averaging over the angles is a well-known approximation technique; in the literature, 
the error of this method has been discussed mainly from a qualitative viewpoint, even 
in the simple case of one frequency (i.e., one angle only). The classical, qualitative 
estimates for this case (see e.g. [5]) have the form 

l{t) - J{et) = 0{e) for t G [0, 0{l/e)) (1.1) 

(uniformly in t) for a perturbation proportional to a parameter e, in the limit e 0^; 
here, l{t) are the actions at time t, and J{et) is their approximation obtained from 
averaging (see paragraph lA for more details). 

In a previous paper [1], we have proposed in place of (11.11) a fully quantitative 
error estimate for the one-frequency averaging; this has the form 

\l{t) - J{et)\ ^ en{et) for t e [0, U/e) , (1.2) 

where | | is the Euclidean norm on the space of the actions and u is a computable 
function, determined by an integral inequality; f/ is a specified nonnegative constant, 
defining quantitatively the time interval where the estimate holds. (In fact, in some 
special cases considered in [1] the estimate holds even for very large values of U, 
e.g., U ~ 1/f^)- Let us repeat here a comment already done in the cited work: 
the idea of a really quantitative approach to the averaging methods has attracted 
little attention up to now, a notable exception being [3] that, in [I], we have briefly 
compared with our approach. 

The present work is an improvement of [1] proposing more detailed error esti- 
mates, e.g., a separate bound on each component of the actions. These componen- 
twise bounds are seen as a special case of a more general framework, where the 
estimates are expressed in terms of any separating family of seminorms on the space 
of the actions (a notion to be defined in the sequel). Our general estimates will take 
the form 

|l(t) - J(£t)|^' ^ £n^(£t) for /i G M, t G [0, U/e) , (1.3) 

where (| |^)^6Af is the family of seminorms, labeled by a (finite) index set M. 

To show the effectiveness of these bounds, in the present work we give a simple 
example related to rigid body dynamics. A more engaging application, concerning 
the motion of a satellite around an oblate planet, will be presented in the companion 
paper [2]. 
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The forthcoming paragraphs lA-lD introduce the following topics: the setting 
of [1] for one-frequency averaging, that we use partly in this paper; the new error 
estimates developed in the present work; the motivations to consider these refine- 
ments, and to formulate them in the language of seminorms; the organization of the 
paper. 

lA. One-frequency averaging, in the framework of [1]. We consider an open 
set A of R*^ and the one- dimensional torus T (referred to as the spaces of the actions 
and of the angular variable): 

A = {/ = {I%=i,...,a} C , T := R/27rZ = {^} . (1.4) 

We suppose to be given a one-frequency system with a perturbation ef on the actions 
and eg on the angle: more precisely, we have a Cauchy problem 



dl/dt = ef{l,Q) , 1(0) = Jo, 

dQ/dt = u;{l) + eg{l,e) , 0(0) = 



:i.5) 



under the assumptions 

/ = {f%=i,...,d e C"^(AxT,R'^), g e C™(AxT,R), cu E C""(A,R) (m ^ 2) , (1.6) 

uj{I) ^ for all / G A , Jq G A, G T , e>0 ; 
the maximal solution (in the future) of (11. 5p is a (7"^+^ function 

(I, 6) : [0, T) ^ A X T , (I(t), e(t)) . (1.7) 

Throughout the paper, the initial data Io,'&o and the perturbation parameter e are 
fixed; for this reason, we do not indicate the dependence of (I, 9) and other functions 
on these objects. Needless to say, we are mainly interested in the case of small e. 
The averaged system associated to (11.51) is 

^=7(J), J(0) = /o, (1.8) 

7 = {T%=i,...4 e c"^(A, R-^) , / ^ /(/) ■■=^Jj^ fii, ^) ; 

the maximal solution (in the future) is a C"^^^ function 

J:[0,W)^A, r^j(r). (1.9) 

The error of the averaging method is the function t ^{t) — J{et) (defined whenever 
l{t) and J{et) exist); equivalently, one can consider the function 

L : t L(t) := i[l(t) - J(rf)] . (1.10) 



In [T] we have put the attention on the Euchdean norm 

d 

Y,u{ty- (1.11) 

i=l 

under natural conditions, we have derived for it a quantitative estimate 

|L(t)| ^ n{£t) for t e [0, U/e) (1.12) 

(which is the same as (11.21) ). where n : [0,U) [0,+oo) is a function determined 
by a fully explicit algorithm. To compute n, one must solve an integral inequality 
or a related differential equation on [0, U), a task that in typical cases is performed 
numerically; however, for e small this operation is much faster than the direct nu- 
merical solution of the perturbed system (11. 5p for t in the long interval [0, U/e). 

IB. Some variants in analyzing L(t). In view of applications, the following 
variants can be of interest: 

(a) estimating a norm of L{t) different from (ll.lip : 

(b) giving separate estimates on the absolute values |L*(t)| of the components; 

(c) considering a partition ^ = {5*, 5"...} of {l,...,d} into (nonempty) subsets 
S,S',... and estimating the components of L(t) in each subset: for example, one 
could analyze the quantities 

y ies y ies' 

Here are some reasons to study each component separately, or to group them into 
subsets: the components could measure physically nonhomogeneous quantities; one 
expects relevant differences in their numerical values, even in the orders of magni- 
tude. All these facts will occur in the example of Section [31 related to rigid body 
dynamics. 

IC. General estimates for L(t) via seminorms. A unified way to treat (a) (b) 
(c) and other situations is to consider on a separating family of seminorms, and 
use them to estimate L(t). Let us recall that a seminorm on R'^ is a map 

R'^^[0,+oo), X^\X\, (1.14) 

homogeneous and subadditive: 

|AX| = |A||X|, + ^ |X| + |r| forX,FGR^AGR (1.15) 

(|A| is the absolute value of A; the first relation, with A = 0, gives |0| = 0). An 
example of a seminorm is the function | |* on R*^, where i is any integer in {1, c?} 
and 

\X\':=\X'\ (1.16) 
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for all X = {X^, ...,X'^) E R*^; more generally, if S* is a (nonempty) subset of R'^ we 
can define a seminorm | on R'^, setting 




A norm on R'^ can be defined as a seminorm with the supplementary separation 
property 

|X|=0 X = 0. (1.18) 

Clearly, this property is lacking (for d > 1) in the example (11.161) : it is also lacking 
in (11.171) . unless S = {1, ...,d}. However, these examples with variable i or S, and 
other situations, carry to families of seminorms possessing the separation property 
in a collective sense. To be precise, a separating family of seminorms on R'^ is a 
family (| |'^)^gM, where M is a finite set, such that | l'^ is a seminorm for each fi and, 
for all X e R^ 

|X|'' = for each /iGM ^ X = 0. (1.19) 

An example of a separating family is formed by all the seminorms (11.161) . with i 
ranging in {1, ...,d}. Another example is the family (I1.17p . labeled by the subsets S 
in a partition of {1, ...,d}. 

Throughout the paper, our estimates for L{t) will concern the nonnegative quan- 
tities 

|L(i)|^ (/iGM) (1.20) 

for any chosen separating family of seminorms on R'^. Case (a) of the previous 
paragraph corresponds to the choice M = {1} and | |^ = a norm | | on R'^; case 
(b) corresponds to the family (I1.16P with M = {1, ...,d}, and case (c) to the family 
(fTTTl) with M = ^. 

ID. Organization of the paper. Section [2] is the main body of the paper: after 
recalling a basic Lemma from [JLJ , we construct the general framework to estimate L 
through a separating family of seminorms. The conclusion is a set of inequalities 

|L(t)|^ ^ n^{et) for E M , t e [0, U/e) (1.21) 

(i.e., of the form (11.31) ). where the estimators : [0, U) [0, +oo) are determined 
solving a system of integral inequalities (Proposition [2]6]), or of differential equations 
related to them (Proposition 12. 7p . Section [3] presents an example, arising from the 
dynamics of a rigid body under damping; this was introduced in [1] and will be 
reconsidered from the present viewpoint, deriving separate error estimates for each 
one of the two actions. The Appendices El [B] contain the proofs of the previously 
mentioned Propositions. 

In spite of the frequent reference to [1], in writing the present paper we have 
tried to make it reasonably self-contained. 
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2 Main results. 

2A. Notations, (i) Throughout the paper, vectors of R'^ are written with upper 
indices: X = ... as already done in the Introduction. Due to the fact that 

R*^ has a canonical basis, the tensors on R"^ of any type (p, q) can be identified with 
tables of real numbers, that we write in the usual style with p upper and q lower 
indices. In the sequel we will often use the tensor spaces 

TliK") = = (i^/) eR for I, J = 1, d} , 

TliR'^) = {^= {^^) I e R for i,j = 1, ...,4 , (2.1) 

T^(R'^) = K = C^y I <^*,GRfor^,j,A; = l,...,4 . 

We use systematically Einstein's summation convention on repeated upper and lower 
indices. Let X,Y e R^ ^ G T}(R°') and ^ G T^(R'^); then, G R"^ and 
G R'^ are the vectors of components (^X)* = i^/X^ {"rfXYY = 'rfjf^X^Y''; 
G T}(R'^) are the tensors of components = (^X)i, = 

%X^ . 

(ii) We fix on R'^ a separating family of seminorms 

I r if^eM) , (2.2) 

with M a finite set. To go on, we need some seminorm families on the tensor spaces 
T}(R^) and T^(R'^); of course a seminorm on TJ(R'^) is a homogeneous, subadditive 
map 

I I : T}(R'^) ^ [0,+oo) , £/^\£/\ (2.3) 

and a seminorm on T2(R'^) is defined similarly. 

Keeping fixed the family (12.21) . a consistent family of seminorms on T}(R'^) is 
one of the form 

11^: (/^,^eM), (2.4) 

with the property 

^ IXI" for all ^ G Ti(R'^), X G R"^ and G M (2.5) 

(here and in the sequel, Einstein's summation convention is also employed for re- 
peated indices with values in M) . Similarly, a consistent family of seminorms on 
T2(R'^) is a family of seminorms 



I UK 



(/i,z/,«:GM) , (2.6) 
such that 

l^XYlf" ^ \'^\^^^\X\''\Y\^ for all G Tj(R'^), X, F G R'^ and /i G M . (2.7) 
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The existence of such consistent famihes can be proved using the separation property 
of (12.21) (0). Of course, if we use on R'^ the seminorms | |* of Eq. fll.l6p we have on 
Tj(R^) and T2(R'^) the following consistent families of seminorms, also taking the 
absolute values of the tensor components: 

|^|}:=K^|, m}fc:=|^,| (z,j,A: = l,...,ci) . (2.8) 

(iii) In the sequel we intend 

At := {(/, (5J) G A X R'^ I [I, I + 5I]cA} , (2.9) 

where [/, / + 61] is the closed segment in R^ with the indicated extremes. 

2B. The integral equation for L. We consider the perturbed and averaged 
systems (ll.Sp (II. 8p . for fixed e > and initial data /q, "do. We introduce the functions 
s G C™(A X T, R^) and p G C'"-i(A x T, R'*) such that 

/^7 + .|i. .^0; (2.10) 



these equations, and the forthcoming ones are always understood in the tensorial 
sense (p). 

From now on, U stands for an element of (0, +oo]. 



2.1 Lemma. Suppose the solution J of ( fi.<§]) exists for r G [0,U). Denote by 
R : [0, U) Tl{R'^), T ^ R(r) and K : [0, U) ^ , t ^ K(r) the solutions of 

dR Of , , , , , , 

- = ^(J)R, R(0) = 1, ; (2.11) 

dK df , ^ .X /X /X 

^ = ^(J)K + P(J) , K(0)=0 (2.12) 



""^This follows from much more general results on multilinear maps and seminorms that can be 
found, e.g., in [5]. 

^For better clarity, let us give only some examples. The equivalents in components of Eq. (|2.10p 

are 

In the forthcoming Eq.s (|2.18p and (|2.20p . the relations about for and /, ^ mean, respec- 

tively: 



Til + SI) = Til) + + l-^kil, 6I)SPSI'^ . 
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(these exist and are C™; R(r) is an invertible matrix for all t G [0, U), and K(r) = 

R(r)/;dr'R(r')-ip(j(r')). Ford=l, R(r) = exp /; dr' |^(j(r')) G (0, +00);. 

Furthermore, assume that the solution (I, 0) of the perturbed system M.5\) exists 
for t G [0, U/e), with {J{et), - J{et)) G Af. Finally, define 

L : [0, U/e) ^ , t^ L{t) := ^ [l{t) - J{et)] . (2.13) 

Then, fort G [0,U/e), 

L{t) = s{l{t), e{t)) - Riet) silo, ^0) - K(£t) (2.14) 

~s(wm,e{t))-^{j{et))v{l{t),e{t)) 

+ e'Riet) ^*rft'R-^(et')(«(I(O,0(O)-|y(J(^O)(«^ + g)(I(O,0(O) 

-^{J{et'))v{l{t'),e{t')) - '^{J{et'),eL{t'))L{t') + ^J^{J{et'), eL{t'))L{t'y 

In the above, v G C"^(Ax T, R'^), q,w e C™"i(Ax T, R'^),ue C^'-^Ax T, R^) and 
M G C™^^(A, r}(i2'^)) are the functions uniquely defined by the following equations: 

dv 

s = uj— , v{I, ^0) = for all I G A; (2.15) 
dw 

P = P + ^1^^ w{I,^o) = for all I eK- (2.17) 

Furthermore, ^ G C™-2(A,, r}(i2^)) an(i ^ G C'"-2(A^, ^^(i^'^)) are two functions 
such that, for all (/, 61) G A|, 

p{I + 61) =p{I)+^{I,SI)SI , (2.19) 

/(/ + SI) = /(/) + ^(/)5/ + -^(/, 5/)5/2 , ^U^, SI) = ^k^I, W) . (2.20) 
Proof. See [1]. □ 
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2.2 Remark. In dimension d = 1, Eqs. fl2.19p fl2.20p can be uniquely solved for 
'^^{1,61) and J^{I,6I); in any dimension we have the solutions given by Taylor's 
formula, i.e., 

dp d^l 
^(1,61) := / dx^(I+x6I), jr(I,6I):=2 / dx (l-x)-^(I+x6I) . (2.21) 
Jo 91 Jq dP 

If p (resp. /) is a polynomial or rational function of the actions, ^ (resp. M') can 
be obtained in a simpler way by direct inspection of Eq. f l2.19p (resp. f l2.20p ). 

Now, from the integral equation fl2.14p for the function t ^— > L(t) we wish to infer 
a system of integral inequalities for the functions t ^— |L(t)|'^, where (| l'^) is any 
separating family of seminorms on R'^. This requires a set of auxiliary functions, 
estimating several characters in fl2.14p . which are introduced hereafter. 

2C. New auxiliary functions. For each set Z, we write 

Z*^ := [z = {z^)^^j^, \z^eZ \^^} . (2.22) 
For J e R'^ and q = (^^) G [0, +cx)]*^ we put 

B{J, q) := {/ G R'^ I |/ - J\^ < V/i G M} . (2.23) 
We further assume the following. 

(i) p = (pt') G ^([0, U), [0, +oo]^0 is a function such that 

fi( J(r), p(r)) C A for r G [0, U) . (2.24) 

We put 

:= {(r, r) G [0, U) x [0, +oo)*^ | r^" < p^(r) V/i G M} . (2.25) 

(ii) a^, b^, c^, d'^, e'^^ G C(Fp, [0, +C)o)) (yU, u, k E M) are functions such that for any 
r G [0, f/), 5J G 5(0, p(r)) and ?9 G T, 

|s( J(r) + 5 J, ^) - R(r)s(/o, ^o) - K(r) ^ ^ a^(r, |5J|) , (2.26) 
\w{^{t) + 5J, d) - |^( J(r)) t;( J(r) + 5J, ^) ['^ ^ 6^(r, |5J|) , (2.27) 

|M(J(r) + 5J, ^) - |y(J(r))(«; + g)(J(r) + 5 J, (2.28) 

-^(J(r))i;(j(r) + 5J,i9)|^ ^ c^{t,\5J\) , 

|^(j(r),5J)|:;^<(r,|5J|), (2.29) 
|jr(j(r),5J)L^ei;,(r,|5J|) . (2.30) 
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In the above, one always intends 

\6J\ := {\6J\^)xeM . (2.31) 

The functions c^id'^i e(^^ are assumed to be nondecreasing with respect to the variable 
r, i.e., 

(r, r), (r, r') G T^, ^ r'^ V A G M ^ c^(r, r) < c^(r, r ) (2.32) 
and similarly for df^ and e(^^. Given a'^, e(^^, we define the functions 

a'^ G C{Tp, [0, +oo)), a^(r, r) := a^(r, r) + £6'^(r, r) , (2.33) 
7^^ G C{T, X [0, +oo)^^ [0, +oo)), (2.34) 

7'^(r, r, £) := c^(r, r) + rfj:(r, T)t + ^e:;,(r, r)rr . 

In the sequel we will set a := (a^) G C{T [0, +oo)*^), and intend 7 similarly, 
(iii) P^^ G C([0, f/), [0, +00)) (/i, z/ G M) are functions such that, for r G [0, f/), 

|R(^)li:^^::(r) , |R-^(r)|:;^P,^(r) . (2.35) 

2.3 Remarks, (a) In the main following statements about |L(t)|'' (Propositions 
12. 6[ 12. 7p . the functions h^, e^^ will always be multiplied by the small factor e. For 
this reason, in applications one can determine 6^, ...,e(^^ via fairly rough majoriza- 
tions of the left-hand sides of Eqs. fl2.27p - fl2.30p . The situation is different for the 
functions a'^, that are not multiplied by e and so require accurate estimates. 

(b) A trivial choice for the functions in (iii) is -R(^(r) := |R(r)|(^, PI}{t) := |R~^(r)|(^. 
This is not satisfactory if one wants more than the regularity: in fact, this choice 
does not grant R'^ and Pj^ to be for any 1. On the other hand, regularity 
with = 1 or 2 will be required by some subsequent manipulations, and in view of 
this we leave R'^ and unspecified. 

2D. Integral inequalities for (jLl'*). We keep the assumptions and notations of 
the previous paragraph. 

2.4 Lemma. Assume that the solution (1,6) of the perturbed system exists on 
[0, U/e) and that |L(t)|^ < p^'{et)/e for all pE M,te [0, U/e), Then, for all and 
t as above, 

\L{t)r^a'^{et,e\L{t)\) + 6'R^,i6t) f dt'P^{6t')r{et',e\L{t%\L{t')\) , (2.36) 

Jo 

intending |L(t)| := {\L(t)\'^)i,<zM ■ 

Proof. We take the /x-th seminorm of both sides in Eq. fl2.14p . To estimate the 
right-hand side, we use the consistency inequalities 02.51) fl2.7p . together with the 
relation | dt'..\^ ^ dt'\..\^] next, we apply the inequalities fl2.26l) - fl2.30l) with 
6J = l(t) — J{et) = eL(t), and the inequalities (12.350 . In this way we obtain 
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\L{t)\^ i^a^'{6t,e\L{t)\ ) + e¥\et,6\L{t)\)+e^R>'^{et) [ dt' P^{et') (2.37) 

x(c^{et',eUt')\) + d:{et\e\m\) UtT + \e:,{et\e\L{t')\) |L(t')r|L(t')r) • 

Now, the thesis fl236|) follows from the definitions ^^,^M) of a, 7. □ 

2E. A general fact on integral inequalities. This result is stated without proof, 
being a simple variation of similar ones appearing in p] [6J. 

2.5 Lemma. Let T e (0,+cx)], S = (S^) E C([0,T), [0,+cx)]^^) and 

E := {{t, e) e [0, T) X [0, +00)^^ I r < S^'it) ^fi}, (2.38) 

E:={it,t',i) \te[o,T), t'e[o,t], {t',i)eE} . 

Consider two funcUons ^ = {^^') e C(H, [0, +00)^'^ andri = [r]'') e C(H, [0, +00)^^). 
Let each function rj^ be nondecreasing in the last variable: ri^{t,t',i') ^ ri^{t,t',i) 
if it,t',£), {t,t',i') e H and ^ t for all u e M; furthermore, let I = (P),D = 
(d^) G C([0,T), [0, +00)^^) be such that graph I, graph C S, and 

no) = , nt) ^ e(t, l{t)) + f dtYit, t', l{t')) , (2.39) 

D'^(t) > e(t, D(t)) + / rftV(i, i', o(^')) (2.40) 
/or a// fi G M, t G [0, T). Then, for all such fi and t, 

['^(t) < o'^(t) . (2.41) 

2F. The main Proposition. We still assume that the solution J of the averaged 
system exists on [0,U), and define R, K via Eqs. (12. lip (12.121) . Moreover, let us 
be given a set of functions p^, a^, 6^, c'^, rf^, e(^^ as in paragraph 2C; and 7^ are 
defined consequently, as indicated therein. 

2.6 Proposition. Assume there is a function n = (n^) G C([0, [/), [0, +00)^^) 
such that, for all fi E M and r G [0, U), 

n^(r) < p^{r)/e , (2.42) 

n'^(r) > a^(r, £n(r)) + £i?^(r) rfr'P,\r') Y{t', 5n(r'), n(rO) . (2.43) 

Then, the solution (1,0) of the perturbed system exists on [0,U/e); furthermore, 
defining L as in Eq. Ii2.13\) we have 

\L{t) 1^ < n^{et) for all e M , t E [0, U/e) . (2.44) 
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Proof. It is given in detail in Appendix \^ however, here we sketch it in few hnes. 
The main idea is to compare the inequahties fl2.36p for |L(t)|'^ and fl2.43p for n'^, 
writing the second one with the change of variables r = et, t' = et' . The thesis 
follows using Lemma [2.51 with := |L(t)|'^, t>^{t) := n'^let) and obvious choices 
for ^'^j i]'^; this is combined with a continuation principle for ODEs, to prove the 
existence of (I, 6) for all t G [0,U/e). □ 

2G. A differential reformulation. We keep the assumptions at the beginning 
of the previous paragraph, but we require some more regularity on the functions 
a^, ...,e(;^, Pl^, Ri^ fulfilling Eqs. (12:261) - fl230D and flOSD . namely 



e C'iT,, R) , e;;, G ^^(r,, R) 



(2.45) 



2.7 Proposition, (i) Assume there are = (^^) G [0 + oo)^'^, (A^) G [0, +oo)^^' 
and a = (a^) G (0 + oo)*^, such that 



A := max A'; < 1/e , 



(0,££) 



^ for e M, i e^, 



\a^{0, eQ - CI + eA'^a" < a'' for fieM. 
Then, there is a unique io = (^q) such that 



(2.46) 
(2.47) 

(2.48) 
(2.49) 

(2.50) 



(ii) With io as above, let m = (m''),n = (n'^) G C^{[0,U), R^^) solve the Cauchy 
problem 

-^ = P>!j^(.,en,n) , m^(0) = , (2.51) 



dn^ / da, f da^ , , „^ ^, , dRl ,^ 



dr 



dr 



X \ dr 



dr 



n^(0) = 



(fi E M) with the domain conditions 

<n'' < pf'/e , det (^l-£^(-,en)) > 



(2.52) 



(2.53) 
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(in the above, 1 — eda/dr stands for the matrix {6^ — eda^/dr^) G M), and 
Eq. \2. 52\) contains the matrix elements of its inverse. Note that Ii2.51\) implies 

^ o;. 

Then, the solution (1,0) of the perturbed system exists on [0,U/e) and 

\L{t)\^ ^ n^(et) for all e M , t e [0, U/e). (2.54) 

Proof. See Appendix [HI also containing a preliminary Lemma. □ 

2.8 Remark. The previous Proposition mentions io, the unique fixed point of the 
map a{0,e-) in the set E of Eq. fl2.46p . The proof in the Appendix indicates that 
a{0,e-) : S — > S has Lipschitz constant eA < 1 in the maximum component norm 
||2;|| := max^ l^;'^!, with A as in fl2.47p . So, from the standard theory of contractions, 
we have the iterative construction 

Eq = lim In, /i any point of S , Z^j := a(0, e/„_i) for n = 2, 3, ... . (2.55) 

n— >+oo 

For each n ^ 2, 

lKo-/„|K(.yi)"-^|^H^. (2.56) 

2H. Implementing the scheme in a typical case: the "!N-operation". Let us 

consider a situation in which (for given data (Jq, ^o) and £ > 0) we have analytical 
expressions for the solution J of the averaged system (11. 8p (on an interval [0, U)) and 
for all the auxiliary functions R,K,s,p, ...,W,J^, a'^, ...,e(^^, Pj^, of the previous 
paragraphs (having chosen a separating family of seminorms (| |^)^gM, and making 
the regularity assumptions (12.451) ). One can provide nontrivial examples where these 
expressions can be obtained: one of them is considered in Section [3l 

In this situation, to obtain the final estimates |L(t)|^ ^ n'^{et) of Proposition 12.71 
we need: the fixed point Iq, defined by Eq. (12.500 : the functions m = (m''),n = 
(n'^) fulfilling the Cauchy problem fl23TD fl232|) . Typically, to find 4 and m'^,n'' 
analytically will be difiicult or impossible, and a numerical approach will be required. 
Concerning £o, one can compute numerically the iterates I2, 13, ...In in f l2.55l) up to a 
sufficiently large order n, and then approximate io with /„. As for m^, n'^, one can 
attack the Cauchy problem (I2.5ip (12.520 by any package for the numerical integration 
of ODEs (paying attention to the domain conditions (I2.53P ). 

From now on, the term "N-operation" will be employed to indicate the numerical 
determination of £0, in, n along the above lines (jf|). Generally, the K-operation to find 
io and m, n on the interval [0, U) is faster (and more reliable) than the computation 
of L(t) := [l{t) — J{et)]/e on the long interval t G [0, U/e), through a direct numerical 

''this is somehow different from the " D^-operation" of [IJ, that also included the numerical 
computation of J, R, K on [0, U) and was designed to work with a single norm on R''. 
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attack to the perturbed system (11.51) : we think that this gives a practical value to 
the general framework developed here. This situation will be exemplified in Section 



21. The "^-operation". This expression will be used to indicate the direct numer- 
ical computation of L from the perturbed system (11.51) . on the time interval [0, U /e); 
in the general framework of this paper, this operation must be performed only if one 
wants to test the efficiency of the Tsf-operation. 

Let us clarify the previous statements, assuming again to have the analytical 
expressions of all the functions mentioned in paragraph 2H. Having the expression 
of J, we substitute I(t) = 3{et) + eL{t) in Eqs. (II. 5p for (I, 6); this gives rise to the 
Cauchy problem 



for the unknown functions t ^ (L(i), ©(i))- By definition, the "XL-operation" is the 
numerical solution of fl2371) for t G [0, U/e) (0). 

The efficiency of the 3\f operation is tested via L comparing: (1) the CPU times 
^L-, 2^?\f required to perform both operations on standard machines; (2) the graphs of 
the functions [L'^l and of their estimators n'^, made available by the two operations. 
Of course, the test is satisfactory if: 

(i) T:k is considerably shorter than T^,; 

(ii) for each e M the estimator 1 1-^ n^(et) approximates well the envelope of the 
rapidly oscillating function t ^ |L(^)|'^, for t G [0, U/e). 

The whole procedures concerning JSf and L are illustrated in the next section; in 
the example therein, both (i) and (ii) will occur in the test of TsT via L. 



{dL/dtm = f{J{et) + eL{t), 0(t)) - f{J{et)), 
{de/dt){t) = cu(J(£t)) + eg{J{et) + eL{t), Q{t)) 



L(0) = 
6(0) = 







(2.57) 



^This differs from the " £-operation" of [T, which included a preliminary numerical determina- 
tion of J on [0, [/). 
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3 An example from rigid body dynamics. 



3A. Introducing the example. We consider a perturbed integrable system of the 
form fll.Sp . with 

d = 2, A := {/= I G (0,+oo)}, w(/) = /^/^ (3.1) 

f{I,^) := (-/i(Ai + /iCOs(2?9)),-/2(A2-/icos(2^)) , g{I,^) :=/isin(2i?) ; 
this depends on three real coefficients /i, Ai, A2 such that 

Ai > , -Xi< fi< Xi , A2 > -Ai . (3.2) 

This system has already appeared in p^j (Section 4, Example 4) where it was related 
to Euler's equation for a rigid body with gyroscopic symmetry under a damping 
moment proportional to e, with a particular dependence on the angular velocity. 
The actions 1^,1^ have different physical meaning: in fact, is the equatorial 
angular velocity (in suitable units) and P measures the inclination of the angular 
velocity on the gyroscopic axis (jfl). We will take for (11.51) the initial conditions 

/o = (/o\/o')eA, ^o-=0. (3.3) 

The forthcoming analysis shows that, depending on the data and on the other pa- 
rameters involved in the problem, the numerical values of the solution components 
I* {i = 1, 2) can be very different over long times; the same happens for the compo- 
nents J* and V{t) = [V{t) — J^{et)]/e. For this reasons, and for the different meaning 
of the two actions, it can be of interest to derive separate estimates for the absolute 
values |L*(t)| {i = 1, 2); this will mark a difference with the analysis of where we 
only gave a global estimate for a/ (L^)'^ + (L^)^. 

3B. Analysis of the example. The average / and the solutions J of (11. Sp . R, K 
of (12. lip (I2.12P (on any interval [0,U)) are written in the forthcoming Table 1, 
which also reports the auxiliary functions s,...,J^ required by our method. As 
anticipated, our aim is to estimate separately the absolute values |L*| (i = 1,2); this 
marks the difference with respect to pQ, where this example was treated estimating 

the Euclidean norm ■\J (L^)^ + (L^)^. In the language of Section [H analyzing the 
components of L corresponds to use the seminorms (11.160 . i.e., 

I P : R2 ^ [0, +00) , X^\X\' := \X'\ {i = 1, 2) . (3.4) 

Whenever necessary, we will use for TliR^), T^(R^) the consistent semmorms l lz.ali . 

The second half of Table 1 contains the functions p\ a\ dj, e'jj^, i?*, Pj k = 
1, 2) required by the general framework of the previous section. The choice of is an 

^See Eq. (4.18) of the cited work. 



14 



obvious consequence of the form of A; a*, 6*, c* have been computed binding the left- 
hand sides of Eqs. fl2.26p - fl2.28l) by elementary means (similar to the ones employed 
for the Euclidean norm estimates of [1], but here applied to the components); we 
observe that a*(r, r) is just a bound on |s(J(r), 5J)|* (for \5J\'' = r'^), because 
K(r) and s(/o,'(?o) in the left-hand side of fl2.26p are zero. The functions d^j,G^jk ci-re 
identically zero due to the vanishing of ^, Jif; the expressions for i?* and are just 
the ones of the matrix elements of R and R~^. 

3C. Results. Starting from the functions in Table 1, the !N"-operation has been 
performed for two choices of the initial data /d,/o ^^d of the parameters U, e, 
li, [2,/i, producing as a main output the estimators n = (n*)i=i,2 on [0, [/); the £- 
operation has been performed as a test, producing directly the function L = (L*)i=i,2- 
Both operations were carried over on a PC, using the MATHEMATICA package. 

The results are summarized in Figures a,b,c,d. These report the CPU times 'Zj^, 
Ti: (in seconds) and allow to compare the actual values of |L*(t)| with our estimators 
n*(r), for t = r/e and r G [0, U). 

It turns out that Tj^r/T^, ~ 1/12 for the first choice of the parameters (Figures 
a, b) and Ij^fZ^ — 1/600 for the second choice (Figures c, d), where the time scale 
U/e is overwhelmingly long. In all cases the functions r 1-^ n^i^) practically coincide 
with the envelopes of the oscillating functions r 1-^ |L*(r/e)|, for r G [0, U). 
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Table 1. A list of functions for the example. 



For I = P) e (0, +oo)^ 1? e T and 57 = (57^ 5P) E +oo) x {-P, +oo) 

sm(2^) f X2 + l^ cos(2^) Ai + 3/x cos(2i^) ^^ 
= ^ ( > 71 p(7) = (0,0), 

usin^^? / 2A2 + 2/i + Ai + u cos(2^?) A2 + 2u + 2Ai + 3u cos(2^?) 



2/1/2 \^ 72 ' 71 

,^ ^, itsin^i? / A2 + «cos^i? Ai + 3jL(C0S^'j?\ ,^ usin^'j? 

«^(^'^) = ^TTTT^ r2 ' n ' «(^'^) 



2/1/2 \^ 72 ' 71 ; ' ' ' ' 47172 

4Ai + 6X2PL + 2A2A1 + nXi + /i(4A2 + 3// + Ai) cos(2i?) + 3//2 cos2(2i?) 

3A2/i + 2A2A1 + IO//A1 + 4A? + 3/i(A2 + 5// + 4Ai) cos(2i?) + 15//2 cos2(2t?) 

7^ 

^(7) = diag(-Ai,-A2), ^(7)=diag(-A2,-A2), ^^(7,W) = 0, ^(7,W)=0 
For T e [0,C/): 

J(t) = (71 e-^^\ 72 ; R(r) = diag(e-^l^ e"^^-) , K(t) = (0, 0) ; 

p^(r) J^(r) (i = l,2). 

For T e [0, [/) and r = (r\ r^) e [0, J\r)) x [0, j2(t)) : 



2(j2(r) -r2) ' - V, w 2(J1(t) - r^) ' 

1 |/i|(4Ai + 4A2 + |//|) ^2. ^._ |/i|(4Ai + 4A2 + 3|/i|) 



(Ji(r) -ri)(j2(r) -r2)2 ' 8(J1(t) - ri)2(j2(T) - r2) 

|/.|(16(Ai + A2)2 + 16|/i|(Ai + A2) + 3|/x|2) 



16(Ji(r) -ri)(j2(r) -r2)2 
|/x|(16(Ai + A2)2 + 16|^|(Ai + A2) + 15|/x|2) 



c^(T,r) := 

'^^^^'''1- 16(Ji(r) -ri)2(j2(r) -r2) 

4(r,r):=0; e},(r,r):=0; 

I^{t) := ; P;{t) := e^^^6} (i, j, A; = 1, 2) . 
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0.2 0.4 0.6 0.8 1 0.2 . 4 0.6 0.8 1 

Figure a. /x = l, Ai = 2, A2 = — 1, -^o=4, Figure b. The same parameters as in 

72^4, £ = 10-2, C/ = l. %y^^Omis, Fig.a. Graphs of n2(T) and |L2(r/£)|. 
%c^0.391s. Graphs of n^(T) and \L^{T/e)\. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

Figure c. /i. = l, Ai = l.l, A2 = — 1, -^o=4, Figure d. The same parameters as in 

/2 = 4, £ = 10-3, U = 3. Tk = 0.047s, Fig.c. Graphs of n2(r) and iL^ir/e)]. 
%c^31.3s. Graphs ofn^(T) and \L\T/e)\. 
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A Appendix. Proof of Proposition 



Let [OjV^/e) (with V G [0,+oo]) be the domain of the maximal solution (1,0) of 
(ll.Sp . and put 

U' ■.= mm{V,U) . (A.l) 

Step 1. One has |L(t)|^ < n^'{et) for all fie M,t e [0, U'/e). To show this, we write 
the integral inequality fl2.43p with t = et, r' = et'] this gives 

ret 

n^iet) > a^iet, en{et)) + e^R'i{et) / dt'P^{et') Yiet', en{et'), n{et')) (A.2) 

Jo 

for te[0, U/e), hence for t e [0, U'/e). 

To go on, we use Lemma 12.41 with the constant U therein replaced by U'; 
this gives Eq. fl2:36|l for t e [0,U'/e). Keeping in mind Eqs. fl236D and fICT) . 
we apply Lemma [23] with T := U'/e, 5^{t) := p^{et)/e, i^{t,t) := a^{et,ei), 
Tj^'it^t'J) := e^ R^^{et)P;^{et')-f^{et',e£,i), l^{t) := |L(t)|^, o^(t) := n^(et) (the re- 
quirement ['^(O) = is fulfilled by construction). Lemma 12.51 gives l^^{t) < t5^(t), 
yielding the thesis. 

Step 2. One has V ^ U, i.e., U' = U (thus (I, 6) exists on [0, U/e) and the inequality 
of Step 1 holds in this interval). Indeed, suppose V <U and put 

K := {{t,I) e[0,V/e]xR'^ \ \I - ]{et)\^' ^en^iet) y jj e M} . (A.3) 

This is a compact subset of R x R"^ and (due to Step 1 and (12.24^ ) graph (I, O) C 
fCxTcRxAxT. But i^' x T is compact: so, by the continuation principle for 
ODEs [7J, the solution (1,0) can be extended to an interval larger than [0,V/e). 
This contradicts our maximality assumption, and concludes the proof. □ 



B Appendix. Proof of Proposition 12. 7L 

We begin with a Lemma, holding under the assumptions at the beginning of para- 
graph 2E. Its proof is very similar to the one given in |T] for Lemma C.l, so it will 
not be reported. 

B.l Lemma. Assume there are functions Ug = (n^) G C([0, t/5), [0, +00)^), la- 
beled by a parameter 6 G (0, 6<^], such that the following holds: 
(i) Us ^U for 5^ 0+; 

(a) for all (5 G (0, 5*], /i G M and r G [0, U5) , 

<ir) < P'{r)/e , (B.l) 
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<(r) = 6 + a^{T,ens{T))+eR^^{T) / dr'T^ir'^r' ,ens{r'),ns{r')) ; (B.2) 
(in) for each fixed r G [0, U), the limit 



n(r) := lim ns{T) 

5~>0+ 



(B.3) 



exists in [0,+oo)*^. 

Then the solution (1,0) of U.5\) exists in [0,U/e) and, defining L as in Eq. 

\L{t)\^ < n^(et) for allfieM,te [0, U/e). (B.4) 

We come to Proposition 12 .71 so, we make the assumptions at the beginning of 
paragraph 2G, strengthened by the smoothness requirements fl2.45p for the functions 

Proof of Proposition 12.71 We introduce the norm || || on R^^, the function 
and its generahzations as {5 G [0, +C)o)) setting 



max z"^ 



cto : S — s> R, 
: S — >■ R , 



i^ao{i) ■.= a{0,ei) , 
i ^ as{i) := ao{i) + 6 



(B.5) 

(B.6) 
(B.7) 



Step 1. For each S ^ 0, as is a contractive map with respect to the norm \\ \\. In 
fact, by ([23HD it is 



^a^' 



dr" 



(0,££) 



for all £ G S . 



(B.8) 



Let £,f G S; the equation = ds{day dr''){{l- sY + stj x(r-f'^) 

implies 

- a^s^n ^ eA>t\t - n , (B.9) 

whence 

\\as{e) - as{(.')\ \ ^ eA\\l~e\ \ < . (B.IO) 

The last two inequahties depend on (12.471) : contractivity of as is proved. 

Step 2. There is 5* > such that, for all 5 E [0, 5*], as sends S into itself. In fact, 

for any 6^0 and £ G S, 

l«5(^) - CI = K(^) + <5 - CI ^ K(^) - <iQ\ + KiQ - CI + ^ 

^ eA^^ir - Q + K(4) -FJ+S^ eA^X + l«o(^*) - CI + ^ , (B-ll) 
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where the second inequahty follows from Eq. fIB.Qp with 6 = 0. To go on, we note 
that fl2.49p implies the existence of a 5* > such that 

|a''(0, et) - CI + eA'^a" + 6, ^ for all /i e M. (B.12) 

For 5 E [0,5,] and £ G S, Eqs. flRTTD . flB:T2D imply < a'^, i.e., as{i) G S. 

Step 3. For all 6 G [0, 5*], the map as has a unique fixed point is = (O £ S, which 
depends continuously on S. Existence and uniqueness of the fixed point follows 
from the Banach theorem on contractions; to prove continuity we note that, for all 
5,5' e [0,(5,], 

Wes-is'W = IMes) - as'MW = ||ao(4) + 5-ao(4') -5'll (B.13) 

^ \\aoii5)-aoiis')\\ + \6-6'\^EA\\£s-is'\\ + \S-6'\ , 
the last inequality depending on fIB.lOP with 6 = 0. This implies 

(B.14) 

so the map 5 i— > £5 is Lipschitz, and a fortiori continuous. 

Step 4. Proving the thesis of (i). This follows from Step 3, with 6 = 0. 

Step 5. Proving the thesis of (ii). For any 6 G [0,5=,,], let is be as in Step 3. First 

of all, let us consider the space of real matrices N = (Nj^) ^j^^i^^m , with the norm 

IliVll := 

sup^6RM,2^o ll^^ll/lkll; ^0^6 that the second inequality (IB.Sp . with 
([23ZD, yields 

e\\^{0,is)\\^eA<l. (B.15) 
By a well-known fact on matrices of the form 1 — A^, this implies 

det(l-£^(0,4)) >0 . (B.16) 
or 

From the standard continuity theorems for the solutions of a parameter-dependent 
Cauchy problem, we know that there is a family {Us,vas,ns)s£{o,s*] with the forth- 
coming properties (a) (b): 

(a) for all 6 G (0, 6,], it is = (m^), = (n^) G C\[0, f/5), R*0; furthermore, 

^ = Pi:r{; ens, ns) , m^(0) = , (B.17) 

^ = (1 - {■,ens)) (-, ens) + sRtP;: ens, ns) + e^m^^ , 

<(0) = (B.18) 



20 



0<<<p^/e, det(l -£—(-, en^)) > (B.19) 

(concerning the domain conditions in the last hne, recall ( 1B.16I) : ^ by ( IB. 171) ). 

(b) One has 

Us ^5^0+ U, ^5^0+ m^(t) ^5^0+ m^(t) for all t G [0, U), (B.20) 

where t/, m'^, are as stated in (ii). 

Let us consider the pair 11x5,1x5 for any 5 G (0,5*]. Then, integrating (IB. 170 . 

<{r)= r dr'Pi:{r')Y{r',ens{T'),ns{T')) for r G [0, (B.21) 
Jo 

Furthermore, from Eq. ( IB.lSp and ( 1B.17P we infer 

0=(1- e—{; ^ns))^^ ~ ( ^ ^"^'^ + ^^^^^ ^"^'^ ""'^ + 

= |:(nS' - a''(-,£nj) - ; (B.22) 

therefore, for r G [0, f/), 

<(r) - a^(r, 51x5(7)) - sR^.ir) m'^sir) = <(0) - a^(0, £1X5(0)) - £i?(;(0) rn^O) 

= £^ - a^(0, ££5) = - «^(4) = 6 (B.23) 

(recall the initial conditions in Eqs. (IRTTD fiRlSl) . Eqs. (iRGl) (iRTl) and Step 3, 
giving e^=a^,{is) = ai;{is)+5). 

From Eqs. dEJl]) we see that 1x5 fulfils Eq. (ICT of Lemma [O Due to 

Eq. (IB. 201) on the limit for 6 —>■ 0^, from the cited Lemma we obtain the thesis. □ 
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